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The Loschmidt echo is a measure of the revival occurring when an imperfect time-reversal procedure is ap- 
plied to a complex quantum system. It allows to quantify the sensitivity of quantum evolution to perturbations. 
An initial quantum state \ipo) evolves during a time t under a Hamiltonian H± reaching the state \ipt}- Aiming 
to recover the initial state |V>o) a new Hamiltonian —Hi is applied between t and It. Perfect recover of |^o) 
would be achieved by choosing Hi to be equal to H\. This is not possible in realistic setups, and there always 
appears a difference between Hi and Hi, leading to a non-perfect recovery of the initial state. The forward 
evolution between t and 2t under the Hamiltonian —Hi is equivalent to a backward evolution from t to under 
£^ ■ Hi, which embodies the notion of time-reversal. The Loschmidt echo studies are focused on the cases where 
the dynamics induced by the Hamiltonians Hi and Hi are non-trivial, or sufficiently complex (like that of a 
' classically chaotic one-particle system or a many-body system). 
• i— i \ 

1 Introduction 

The concept of time-reversal has captured the imagination of physicists for centuries, leading to numerous vivid 
discussions. An emblematic example of these was the controversy around the second law of thermodynamics 
between Ludwig Boltzmann and Joseph Loschmidt. When Boltzmann was trying to develop the microscopic 
£*T) ' theory of the second law of thermodynamics, Loschmidt raised an objection that had profound influence on 
• the subsequent development of the theory. He argued that, due to the time-reversal invariance of classical 
' mechanics, evolution in which the entropy can decrease must exist. These states could be reached by reversing 
. velocities of all molecules of the system. After such a reversal the entropy would then no longer grow but 
decrease, seemingly violating the second law of thermodynamics. As a response Boltzmann argued that such 
a time-reversal experiment would be impossible and put forward a statistical interpretation of the second 
^ ' law. Undoubtedly, Boltzmann's argument is a valid approach to the problem of the "arrow of time" for generic 
macroscopic systems. However, in quantum systems with few degrees of freedom, today's technological advances 
make it meaningful to address time-reversal experiments. 



1.1 Pioneering experiments 

Refs. [Hahn 501 iRhim 711 TBrew 841 IZhan 921 UTevs 98] 

The first controlled experimental implementation of time-reversal was achieved in the fifties with the inversion 
of nuclear spins processing around a magnetic field. Such a reversal was able to compensate the local-field 
inhomogeneities responsible for the decay of the free induction signal. The time-reversal was implemented 
through a radio-frequency pulse, leading to the formation of a revival in the induction signal, known as spin 
echo or Hahn echo. Early on Erwin Hahn recognized that his procedure, which he viewed as a change in the 
sign of the system Hamiltonian, provided a quantum implementation of the Loschmidt proposal. Interactions 
between spins, not reversed in Hahn's procedure, were the main cause of the decay of the echo, within a time 
scale known as Ti. This experiment was followed by a number of variants, both within magnetic resonance as 
well in other time-dependent spectroscopies. In particular, the dynamical decoupling techniques implemented 
in quantum registries to isolate them from their environment are variants of Hahn's original experiment. 

The next level of complexity was developed with the reversal of the many-body interactions. No general 
recipe is available in this case, but the sign of the truncated spin-spin dipolar interaction, which is associated 
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with the cosine of the angle between the inter-spin vector and the quantizing magnetic field, can be reversed by 
rotating the spins into a quantization axis. This was the proposal of the Magic Echo procedure, implemented by 
Won-Kyu Rhim, Alex Pines and John Stewart Waugh in the seventies. Substantially simpler is the Polarization 
Echo sequence, implemented by Richard Ernst and collaborators in the nineties. In this last case the initial 
state has a local nature as it is labeled by the presence of a rare 13 C which plays the role of a local probe to 
inject and later on detect the polarization of a nearby 1 H immersed in a *H network. 

This idea was further exploited by Patricia Levstein, Horacio Pastawski and Gonzalo Usaj to test the 
stability of many-body dynamics. They suggested that the inefficiency of the time-reversal procedure found in 
all the previous experiments has a connection with quantum chaos and the inherent dynamical complexity of 
the many-body spin system. 

1.2 Definition 
Refs. jPere 84L IJala 01] 

The Loschmidt echo is defined as 

M(t) = \(i>o\e lH2t/h e~ lH ^ h \tPo)\ 2 (1) 

where 

• \ipo) is the state of the system at time 

• Hi is the Hamiltonian governing the forward evolution 

• Hi is the Hamiltonian governing the backward evolution 

• t is the instant at which the reversal takes place 





Figure 1: Schematic flow of the time-evolution for (a) the Loschmidt echo and (b) the fidelity. 

The time evolution appearing in Eq. ([1} is schematically represented in Fig. HJa) , where the Loschmidt echo 
quantifies the degree of irreversibility. Alternatively, Eq. ([lj can be interpreted as the overlap at time t of two 
states evolved from |i/>o) under the action of the Hamiltonian operators H\ and Hi. In this case the Loschmidt 
echo is a measure of the sensitivity of quantum evolution to perturbations. The quantity M(t) interpreted in 
this manner is illustrated in Fig. [TJb) and is usually referred to as fidelity. The equivalence between Loschmidt 
echo and fidelity is displayed in Fig. [5] for the example of an initially localized wave-packet in a Lorentz gas. 
The probability density of the evolved state under Hi (Hi) is plotted in Fig.[2](b) (c), while Fig. [2jd) presents 
the state resulting from the combined evolution Hi and then — H2 both for a time t. The overlap squared 
between the states (a) and (d) defines the Loschmidt echo, while the overlap squared between the states (b) 
and (c) defines the fidelity. It is important to remark that even though the probability density distributions 
in figures (b) and (c) are seemingly identical, the phase randomization due to the difference between the two 
Hamiltonians leads to a weak Loschmidt echo. 

The definition given by Eq. ([TJ assumes that the Hamiltonian operators Hi and Hi are independent of time. 
A generalization to the case of time dependent Hamiltonians is straightforward. The case of non-Hermitian 
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Gxp[-\H 2 t] ^ exp[i// 2 t] ^ 




Figure 2: Wave-packet evolution in a Lorentz gas. (a) Initial state at time t — with momentum pointing to the 
left, (b) State evolved with the Hamiltonian Hi in the interval (0,t). (c) State evolved with the Hamiltonian 
H2 in the same time interval, (d) State evolved from that depicted in panel (b) with the —Hi for the time 
interval (t,2t). In panels (b) and (c) the classical trajectories corresponding to three initial positions within 
the original wave-packet are shown for reference. The square of the overlap between the states (a) and (d), 
the Loschmidt Echo, is M(t) = 0.09, the same as that between the states of panels (b) and (c). Adapted from 
Ref. |Cucc 02b] . Copyright (2002), American Physical Society. 



operators Hi and H2 can also be considered, but in such a case the equivalence between the Loschmidt echo 
and fidelity would not hold. 



1.3 Loschmidt echo and decoherence 

Refs. jZure OIL ICucc 03L IZure 03L ICucc 04L ICasa 10L IBonaTT] 

In an isolated quantum system the evolution is unitary; initially pure states remain pure. When the system is 
connected to the external world the purity of an initial state is generically lost in the course of time evolution. 
This ubiquitous process is referred to as decoherence. The coupling to the environment degrees of freedom, or 
alternatively the lack of complete knowledge of the system Hamiltonian, typically wash the specific interference 
properties that would be observed within a unitary evolution. Since interference is one of the most important 
signatures of quantum mechanics, decoherence has both, fundamental and practical interest. Concerning the 
former aspect, decoherence has been proposed as a road towards the classical behavior observed in macroscopic 
systems. On the second aspect, it is clear that decoherence represents a limitation in the implementation of 
quantum computers or nanoscopic devices based on quantum effects, specially when the number of qubits is 
scaled up, increasing the complexity of the system. The Loschmidt echo, by including a non-controlled part of the 
Hamiltonian, or the noisy effects of the environmental degrees of freedom, gives a way to quantify decoherence 
effect. By attempting the time- reversal of the controlled part of the Hamiltonian, one can quantify at which 
rate the neighborhood of the system (either a unitary part of larger system Hamiltonian or though the effect 
of infinite uncontrolled environmental degrees of freedom) acts as decoherence. Moreover, other quantifiers of 
decoherence, such as for instance the purity, can be shown, under certain assumptions, to be characterized by 
the same dependence on the underlying classical dynamics as the Loschmidt echo. 



1.4 Broad interest of the Loschmidt echo 

Since the beginning of this century, the Loschmidt echo has been a subject of intensive studies by researchers 
from different scientific communities. A list of problems, in which the concept of the Loschmidt echo appears 
naturally, includes 

• Quantum chaos, or the quantum theory of classically chaotic systems 

• Decoherence or the emergence of "classical world" (open quantum systems) 

• Quantum computation and quantum information 
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• Spin echo in Nuclear Magnetic Resonance 

• Linear waves (Elastic waves, microwaves, time-reversal mirrors, etc.) 

• Nonlinear waves (Bose-Einstein condensates, Loschmidt cooling, etc.) 

• Statistical mechanics of small systems 

• Quantum chemistry and molecular dynamics 

• Quantum phase transitions and quantum criticality 

• Mathematical aspects of the Loschmidt echo 

Such broad interest stems from the fact that the Loschmidt echo is a measurable quantity exhibiting, in 
certain regimes, a robust and universal behavior. It is by means of the forward and backward time-evolutions 
that the Loschmidt echo filters out some uninteresting effects, while amplifying other, more important physical 
processes taking place in complex quantum systems. 

2 Loschmidt echo: techniques and main results 

The Loschmidt echo is typically a decreasing function of t. It is therefore of foremost importance to determine 
the form of the decay and the associated characteristic times in various physical situations as a function of the 
characteristic quantities of the problem. Even though the Loschmidt echo has been addressed in a variety of 
scenarios, the main effort has been directed towards the simplest set-up of one-body dynamics. Well established 
results for the case of one-body systems, together with techniques adapted from the field of quantum chaos, are 
summarized in this section. 

2.1 Characteristic quantities 

A number of physical quantities, describing characteristic time and energy scales of a system under consideration, 
prove to be especially important in the theoretical studies of the Loschmidt echo. These quantities include 

• g, mean density of states. - By definition, g(E) = ^n) d d d^T" 1 ' wnerc 1S the volume of the phase- 
space enclosed by a surface of constant energy E. The mean level spacing A is given by an inverse of 
the density of states, A = 1/g. Note that g is the smooth part (or average) of the full density of states, 
gtu\\{E) = ~ E n ), where {E n } are the eigenenergies of the Hamiltonian. 

n 

• ffiocaii local density of states. - Let be a reference quantum state. The local density of states with 
respect to \ip) is defined as g\ OC a\(E, — \ ( n \i>)\ ~ En), where {|n)} and {E n } are respectively 

n 

the eigenstates and eigenenergies of the Hamiltonian in question. 

• N, effective size of the Hilbert space for a quantum state. - An initial quantum state |?/>o) of a closed 
quantum system can be viewed as a linear combination of N eigenstates of the Hamiltonian. In a d- 
dimensional system, N is given by the volume of the phase space, accessible to \ipo) in the course of its 
time evolution, divided by the volume of the Planck's cell, (27rft) . 

• A, mean Lyapunov exponent. - In classical systems with chaotic dynamics, a distance between two typical 
trajectories, initiating at infinitesimally close points in phase-space, grows exponentially with time. The 
rate of this growth, averaged over all trajectory pairs at a given energy, is the mean Lyapunov exponent 
A. The latter characterizes dynamical instability of a classical system with respect to perturbations of its 
initial conditions. 

• is, Ehrenfest time. - The center of an initially localized wave packet stays close to the phase-space 
trajectory of the corresponding classical particle for times shorter than the so-called Ehrenfest time, t^. 
The latter, in chaotic systems, can be estimated as t-& = j- In — , where A is the mean Lyapunov exponent 
of the classical system, L is the characteristic linear size of the system, and a is the initial dispersion of 
the wave packet. 

• tii, Heisenberg time. - The discreteness of the energy spectrum of a closed quantum system becomes 
dynamically important at times comparable to (and longer than) the so-called Heisenberg time, tn = 
ft/A = hg, where A and g are the mean energy level spacing and the mean density of states, respectively. 
Semiclassical (short-wavelength) approximations to quantum dynamics are known to break down beyond 
the Heisenberg time. 
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2.2 Calculational techniques 
2.2.1 Semiclassics 

Refs. jJala OIL IVani 03L ICucc 04L lyanToil IVani 06l IGutk iol IZamb 11] 

In this context, semiclassics stands for the short-wavelength approximation to the quantum evolution in terms of 
classical trajectories. The key ingredient is the Van Vleck-Gutzwiller approximation to the propagator (matrix 
element of the evolution operator in the position representation) 

v 7 7(q^q',t) v 7 



7 stands for the classical trajectories going from q to q' in time t. 



• Rj — f drC 7 is the Hamilton's principal function for the trajectory 7. Here, £ 7 denotes the Lagrangian 
along 7. 

• C 7 = I det(— V q 'V q i? 7 )| is the stability factor of 7. 

• is the number of conjugate points along 7. 

• d is the number of dimensions of the position space. 

A combination of Eqs. (HJ and (HJ leads to the semiclassical approximation to the Loschmidt echo, 

Oil 



m) = J IK £ liil >*> 

3=1 «i, a 2 , Pi, Pi 



Pi 




(3) 



Here, the summand is a product of terms with the following pictorial representation: 

• A red circles at a point q^ corresponds to (qj\ip). 

• A blue circles at a point q^ correspond to (ip\qj). That is, every blue circle is a complex conjugate of the 
corresponding red circle. 

• A solid curve leading from q., to q/j along a trajectory 7 corresponds to (2-Kih)~ d / 2 C-J/ 2 exp (iR^/h — 
iiris^/2). The Hamiltonian Hi is used if 7 equals ai, and the Hamiltonian H2 is used if 7 equals ^2- 

• A dashed curve leading from q^ to q^ along a trajectory 7 corresponds to (— 27rift)~ d / 2 C^ 2 exp ( — iR~Jh+ 
i7w 7 /2). The Hamiltonian H\ is used if 7 equals 02, and the Hamiltonian H2 is used if 7 equals P\. In 
other words, every dashed curve is a complex conjugate of the corresponding solid curve. 

In Loschmidt echo studies one is typically in interested in the case of H2 being close to H\. It is common 
to describe the perturbation by an operator defined as 

k£ = H 2 -Hi, (4) 

where k parametrizes its strength. 

The expression in Eq. is, in general, extremely complicated to calculate as it involves six spacial integrals 
and the sum over four classical trajectories. Semiclassical evaluations of the Loschmidt echo deal with different 
approximations to Eq. ([3]). The following steps are usually implemented in the semiclassical calculations in 
order to obtain meaningful analytical results: 

(i) If the initial state is spatially localized around a point qo, the diagram of Eq. ([3]) can be approximated 
by a simpler diagram in which the points qi, q2, q3, and q4 merge into qo- Within this approximation 
there are only two spacial integrals to be done (over qs, and q@). 
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(ii) Assuming that the perturbation kS is classically small but quantum mechanically significant (which means 
that the perturbation does not change the topology of the trajectories but introduces a phase difference) 
trajectories a? and j3j are identified for j = 1,2. This is the so-called diagonal approximation. The 
shadowing theorem ensures that the identification of classical trajectories a and j3 is always possible and 
that the resulting action difference can be obtained from the accumulated phase of the perturbation along 
one of the trajectories. For instance, in the case where £ only depends on position coordinates 

Ai? 7 = -k [ dr S(q 7 (r)) , (5) 
Jo 

for 7 = OL\,ot2- The diagonal approximation reduces Eq. ([3]) to a sum over only two trajectories, Oi and 
a2 , of the unperturbed system. 

(iii) A key step that allows to make further progress in the semiclassical calculation is to regroup the resulting 
pair of trajectories (0:1,0:2) into two different families: 

— Trajectories where q5 is near qg and oi evolves close to 02- 

— The rest of pairs, where the two trajectories are uncorrelatcd. 

The semiclassical calculation usually proceeds by estimating the accumulated phases along the trajectories 
for the two kinds of resulting pairs, and then performing some kind of average of M(t) (over the perturbation, 
the initial conditions or the evolution of the classical trajectories). In the case of complex classical dynamics 
these averages are justified and allow to obtain the mean value M(t) of the Loschmidt echo in terms of the main 
parameters of problem. 

An alternative semiclassical formulation for fidelity amplitude m(t) (defined throught M(t) — \m(t)\ 2 ) which 
avoids the usual trajectory-search problem of the standard semiclassics, is the so called dephasing representation 



m(t)= J dqdp cxp (-~AR(q, p, t)j W (q,p 



) (6) 



where 



W o(^P)=(CT yd(«q)exp(-ip.$q) (q + ^|Vo)(^o|q - ^) (7) 
is the Wigner function of the initial state \4>o) and 



AR(q, p, t) = [ dr[£ 2 (q(r), p(r)) - A(q(r), p(r))] 
Jo 



(8) 



is the action difference evaluated along the phase-space trajectory (q(r),p(r)) evolved from (q, p) under the 
unperturbed Hamiltonian H\. In this way the decay can be attributed to the dephasing produced by the 
perturbation of the actions -thus the name dephasing representation. In the case where £ only depends on 
position coordinates (as in Eq. fl5J) the phase difference © reads 

AJ*(q,p,t) = -K / drS(q(r)) . (9) 



Jo 

For generic chaotic systems and initially localized states, the calculation of M(t) using the dephasing repre- 
sentation follows the same lines of the standard semiclassics. The standard semiclassical approach becomes an 
initial-value problem once the integration over the final position is traded by the integration over the initial 
momentum using the stability factor C 7 as the Jacobian of the transformation. However, Eq. (|6]) bears the 
advantage that it regards m(t) as the solution of an initial-value problem, as opposed to a boundary-value 
problem. This proves especially convenient in situations when the sums over classical trajectories are evaluated 
explicitly. 



2.2.2 Random matrix theory 

Refs. jCucc 02aUCerr 031 iGori 041 IStoc 041 IStoc~05l IKohl"08l IKohl HUKohl 12] 

Random matrix theory is a powerful technique to understand the statistical behavior of quantum complex 
systems. Among the latter quantum systems exhibiting chaotic dynamics in the classical limit are particularly 
important, as their spectral properties are well described by averages taken over appropriate ensembles of 
matrices that satisfy certain general symmetry restrictions. The simplest Hamiltonian ensemble is that of real 
symmetric N x N matrices (appropriate for cases with time-reversal symmetries) usually referred to as Gaussian 
Orthogonal Ensemble. 
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The invariance of the probability distribution under orthogonal transformations and the statistical indepen- 
dence of the matrix elements lead to a joint probability distribution of the matrix elements of the Hamiltonian, 



where v defines the scale of the matrix elements and Km is a normalization constant. Each of the N{N + l)/2 
independent matrix element i3y is a zero-centered Gaussian such that 



The averages are taken over the matrix ensemble. The Loschmidt echo in complex systems can be addressed 
by imposing the Random Matrix hypothesis for the forward evolution Hamiltonian Hi, for the perturbation 
kS representing the imperfection of the time-reversal, or for both. Generically, the same results are obtained, 
independently of the choice of the matrix considered to be random. Averages of the Loschmidt echo over matrix 
ensembles are justified for ergodic classical dynamics. The absence of finite classical time scales, like the inverse 
Lyapunov exponent or the escape rate, in the Random matrix theory hinders the description of the Loschmidt 
echo decays that depend on those quantities. 

Recently, using the Random-Matrix- Theory approach, Kohler and coworkers have established relations be- 
tween the averaged fidelity decay and the so-called cross-form factor, characterizing parametric correlations of 
the energy spectra of the unperturbed and perturbed systems. Notably, these relations exist not only in the 
case of fully chaotic unperturbed and perturbed systems, but also in the case of a regular system perturbed by 
a chaotic perturbation. 

2.2.3 Numerical simulations 
Refs. |Tal 84L iDe R 96L rVarTTl] 

A numerical simulation of the time-evolution of initially localized wave-packets under slightly different Hamil- 
tonians, Hi and H2, is a valuable tool in understanding the behavior of the Loschmidt echo in different systems 
and various regimes. Even though there exist numerous approaches to solving the time-dependent Schrodinger 
equation on a computer, the majority of numerical studies of the Loschmidt echo have been concerned with 
the following two methods. Both methods provide accurate, efficient, and stable approximations to the evo- 
lution operator e~ lAr , where the operator A corresponds a (properly rescaled) Hamiltonian, and r denotes a 
sufficiently short propagation time-step. 

• Trotter- Suzuki algorithm. - This method involves three implementation stages. First, one decomposes A 
into a finite (and practically small) number of components, A = A\ + A2 + . . . + A n , such that the operator 
e — iAjt can kg constructed analytically for all j = 1, 2, . . . , n. Second, one defines the symmetric operator 
U 2 (t) = e ~iA n r/2 _ e -iA 2 r/2 e - l A 1 r e - l A 2 r/2 _ e -iA n r/2^ Thirdj one constructs the operator J7 4 (t) = 
U2(pt)U2(pt)U2((1 — Ap)r)U2{pr)U2{jpr) with p = 1/(4 — 4 1 / 3 ). This provides a unitary approximation 
to the original propagator accurate up to the 4th order in r, i.e., e~ lAr = U±{t) + 0(t 5 ). 

• Chebyshev-polynomial expansion. - This approach requires that the operator A is normalized in such a way 
that all its eigenvalues lie in the interval between -1 and 1. The method involves constructing the operator 
Sn(t) = Jo( T ) + 2X)„ = i( _ «)" ' Jn{T~)T n {A), where {J n } and {T n } are, respectively, the Bessel functions 
of the first kind and the Chebyshev polynomials of the first kind. In actual computations, one makes use 
of the recurrence relation T n+1 (A) = 2AT n (A) - T n _i(A), with T (A) = 1 and Ti(A) = A, to calculate 
Sn(t). Then, for large values of n and fixed r, the Bessel function J n (r) rapidly decays as a function of 
n. In fact, |J n (r)| ~ (r/2) n /n! as n — > 00. This is why Sn(t), with a sufficiently large N, provides an 
extremely accurate approximation to the propagator in question, i.e., e~ lAr = Sn{t) + C((r/2) iv /iV!) . 

2.3 Decay of the Loschmidt echo: different regimes 

The decay of the Loschmidt echo mainly depends upon the underlying classical dynamics of the system, the 
initial state, and the nature and strength of the perturbation kE. 

The behavior of the Loschmidt echo is fairly well understood for single-particle quantum systems whose 
dynamics is fully chaotic in the classical limit. Some progress has been done towards the theory of the Loschmidt 
echo in systems with regular and mixed phase space. 




(10) 




i<j, 
(1 + <%)v 2 i<j. 
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Figure 3: Decay regimes of the Loschmidt echo in quantum systems with chaotic classical limit. Characteristic 
crossover times are t c ~ Ji/r], signaling the end of the initial parabolic decay, and t s (~ r _1 TniV for the 
non-perturbative/exponential regime) indicating the onset of the saturation. 



2.3.1 Chaotic dynamics 

Refs. [Jala OIL [Jacq 01[ ICerr~02l IWisn~02l ICerr 03l IWisn 03L ICucc 04L iGori 04L iGuti 09| 

The time-decay of the Loschmidt echo, averaged over an ensemble of initial states, Hamilton operators, or 
perturbations, typically exhibits three consecutive stages (see Fig. 

• Short-time parabolic decay, M(t) ~ 1 — (r/t/h) 2 . It is a short, initial stage of the Loschmidt echo 
decay in all quantum systems. Here, r\ is an average dispersion of the perturbation operator evaluated 
with respect to the initial state, i.e., (t]/k) 2 — (?/>|£ 2 1^) — (V>|£|V') 2 - The parabolic decay holds for times t 
short enough for the propagators of the unperturbed and perturbed systems, exp (—iHjt/K) with j = 1, 2, 
to be reliably approximated by their second order Taylor expansions, 1 — iHjt/H — (Hjt) 2 / (2h 2 ). 

• Intermediate-time asymptotic decay. The short-time parabolic decay is typically followed by an 
"asymptotic" decay regime, whose functional form depends on the strength of the perturbation. 

— Perturbative/Gaussian regime, M(t) ~ exp(—(rjt/H) 2 ). This decay holds for "weak" Hamilto- 
nian perturbations, such that the absolute value of a characteristic matrix element of the perturbation 
operator is small compared to the mean energy level-spacing of the unperturbed Hamiltonian. Since 
77 ~ k, the decay rate is quadratic in the perturbation strength. 

— Non-perturbative/Exponential regime, M(t) ~ exp(— Tt). This non-perturbative decay regime 
is typically observed for stronger Hamiltonian perturbations, i.e., for perturbations large on the scale 
of the mean level spacing of the unperturbed Hamiltonian. The functional form of the dependence 
of the decay rate T on the perturbation strength k is different for "global" and "local" Hamiltonian 
perturbations (see below). 

• Long-time saturation, M[t) ~ TV -1 . At long times, the asymptotic decay of the Loschmidt echo is 
followed by a saturation (or freeze) at a value inversely proportional to the size N of the effective Hilbcrt 
space of the system. N is given by the volume of the phase space, that is available to the state of the 
system in the course of its time evolution, divided by the volume of the Planck's cell, (2Trh) d . The value, 
at which the Loschmidt echo saturates, is independent of the perturbation strength. 

An explicit expression is available for the saturation value of the Loschmidt echo in two-dimensional 
quantum billiards, with strongly chaotic dynamics in the classical limit. Thus, if the initial state of the 
particle is given by the Gaussian wave function ip(q) = (tkj 2 )^ 1 ^ 2 exp (ipo(q — qo)M — (Q — c to) 2 /(2o' 2 )) 
with qo, po, and a denoting respectively the average position, momentum, and position uncertainty of 
the particle, and the area of the billiard is A 7 then the long-time saturation of the average Loschmidt 
echo is given by M(t) ~ ^/2tt ha j '(|po|A) . This formula holds in the semiclassical regime, such that 
VI Pol < o < VA 

The above classification is schematically illustrated in Fig. |U 
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Figure 4: Typical regimes of the Loschmidt echo decay in chaotic systems. 
Global perturbations 

Refs. jJala OIL [Jacq 01| , [C ucc 02bl ICucc "04l IPros 02b] 

A Hamiltonian perturbation is said to be "global" if it affects all (or a dominant part) of the phase-space that 
is accessible to the system in the course of its time-evolution. The accessible phase-space depends on the initial 
state of the system. 

r I 



A 




Kiq ^lo ^ 



Figure 5: Exponential decay rate T as a function of the perturbation strength k in the case of a global Hamil- 
tonian perturbation. 

For global perturbations, the dependence of the decay rate T on the perturbation strength n can be ap- 
proximated by a piecewise-continuous function, shown schematically in Fig. [5] More precisely, the semiclassical 
analysis of the exponential decay regime of the Loschmidt yields 

M(i) = e- CK2t +b(K,t) e~ xt . (11) 

The first term in Eq. (jlip stems from uncorrelated pairs (pti, o^) of contributing trajectories, while the for second 
term is a contribution of correlated trajectory pairs, such that a.\ evolves close to «2 (see Sec. 12. 2.1]) . Unlike the 
constant c, the prefactor b generally exhibits an explicit dependence on both t and k. However, this dependence 
is sub-exponential, and, as a consequence, the decay rate T of the Loschmidt echo is dominated by the minimum 
of ck 2 and A, leading to a commonly accepted (but oversimplified) interpretation of the exponential decay (see 
Fig. IS]). 

For perturbation strengths k that are weak compared to a critical strength k c , the dependence is parabolic, 
r = ck 2 . The semiclassical calculations give the expression of c in terms of correlation functions of the pertur- 
bation, and the random matrix approach leads to the general dependence ck 2 = 2irrj 2 /(HA), where r\ is given by 
the mean value of the width of the local density of states with respect to the initial state IV'o)- Such a behavior 
is commonly referred to as the Fermi-golden-rule regime. 

For perturbations that are stronger than k c , but weaker than a certain bound value Kb, the decay rate T is 
independent of k and equals the average Lyapunov exponent A of the underlying classical system. This decay 
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is known as the perturbation-independent or Lyapunov regime. This regime holds for perturbations 
up to the breakdown strength Kb, beyond which it is no longer possible to regard every trajectory of the 
perturbed system as a result of a continuous deformation of the corresponding unperturbed trajectory, so that 
bifurcations have to be taken into account. The Lyapunov regime is especially remarkable, as the decay rate is 
totally independent of the strength of the perturbation which is at the origin of the decay. 

The exponential decay of the Loschmidt echo is generally followed in time by a saturation at a value on 
the order of TV -1 , where N is the size of the effective Hilbert space. This implies that, in the non-perturbative 
regime, the saturation occurs at time t s ~ In N. In the case of a two-dimensional chaotic billiard (introduced 
above), the saturation time is given by t s ~ r~ x In (jpo\A/(ha)) — (2r) _1 ln(2-7r). It is interesting to note that 
the saturation time can, in principle, be arbitrarily long. In particular, t s can exceed other important time 
scales of quantum dynamics, such as the Heisenberg time tn = %(-£') with g(E) denoting the density of states. 
Indeed, in two-dimensional billiards the saturation time t s — > oo as |po| — > oo, while the Heisenberg time in is 
independent of the particle's momentum. 

The exponential decay, M(i) ~ exp(— Ft), only holds for perturbations weaker than Kb- Beyond this thresh- 
old, the exponential regime breaks down and gives way to another regime, in which the Loschmidt echo exhibits 
a Gaussian dependence on time. 

Some of the mentioned regimes can be obtained using alternative approaches relating the Loschmidt echo 
with the two-point time auto-correlation function of the generator of the perturbation. 

Local perturbations 

Refs. |Gous 071 TGous 081 iHohm 081 [Ares 091 [Kobe llj 

A Hamiltonian perturbation is said to be "local" if it is concentrated in a small region of the phase space 
accessible to the system. In the case of chaotic dynamics, the phase space extent of a local perturbation can be 
characterized by a rate 7, known as the "escape" rate, that is defined as the rate at which trajectories of the 
corresponding classical system visit the perturbation region. For local perturbations, the escape rate is small 
compared to the characteristic rate at which a typical trajectory "explores" the dynamically available phase 
space, and also small compared to the average Lyapunov exponent of the system. 

r a 

a 




Figure 6: Exponential decay rate T as a function of the perturbation strength k in the case of a local Hamiltonian 
perturbation. 

Figure [5] illustrates the characteristic dependence of the decay rate T on the perturbation strength k. The 
function T(k) exhibits a number of distinctive features. Thus, for sufficiently weak perturbations, the decay rate 
grows quadratically with the perturbations strength, r ~ k 2 , demonstrating the Fermi-golden-rule regime. 
In the limit of strong perturbations, k — > 00, the decay rate T saturates at a perturbation-independent value 
27. This saturation is known as the escape-rate regime. The crossover from the Fermi-golden-rule to the 
escape rate regime is non-monotonic, and T(k) generally exhibits well-pronounced oscillations. The amplitude 
and frequency of these oscillations depend on the nature and physical properties of the particular Hamiltonian 
perturbation. An important point is that the function T(k) is given for all strengths k by the width of the local 
density of states. 

The exponential decay, M(t) ~ exp(— Tt) with the non-monotonic function T(k) described above and illus- 
trated in Fig. [6j has been obtained by semiclassical analysis. (Its existence has now been observed in numerical 
and laboratory experiments). The key assumption of any semiclassical reasoning is that the de Broglie wave- 
length corresponds to the shortest length scale of the system. This assumption is obviously violated for systems, 
in which the spatial linear size of the Hamiltonian perturbation is small compared to the de Broglie wavelength. 
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Generally, the decay of the Loschmidt echo in such systems is different from the exponential decay. For instance, 
two-dimensional chaotic systems in the limit of point-like Hamiltonian perturbations, for which 7 — > 0, exhibit 
the algebraic, inverse-quadratic decay of the Loschmidt echo, M(t) ~ i~ 2 . 

2.3.2 Regular and mixed dynamics 

Refs. [Pros 02al iPros 031 [Jacq 03[ [Sank 031 IWem~05l |Wang~07] IGous 11| 

Time dependence of the Loschmidt echo in quantum system with regular or mixed classical dynamics is typically 
more complex than that in fully chaotic systems. As a result of this complexity, a comprehensive classification 
of decay regimes is still lacking. Nevertheless, a number of important results have been established for the 
Loschmidt echo decay in systems whose phase space is predominantly regular in the classical limit. 

In the limit of weak Hamiltonian perturbations, k — > 0, the average Loschmidt echo generally exhibits 
a Gaussian decay. The duration tQ of this decay is inversely proportional to the perturbation strength, 
tQ ~ The existence of the Gaussian decay requires the perturbation to be sufficiently weak such that the 
decay time tQ is long in comparison with any relaxation (averaging) time scale of the system. 

If the Hamiltonian perturbation is sufficiently strong and varies rapidly (or is almost "random") along a 
typical classical trajectory of the unperturbed Hamiltonian, then the average Loschmidt echo is known to exhibit 
the algebraic decay M(i) ~ t~ 3d / 2 , with d being the dimensionality of the system. This power-law decay is 
faster than the decay of the overlap of the corresponding classical phase-space densities, t~ d . 

Numerous case studies of the unaveraged (individual- realization) Loschmidt echo M(t) have revealed that 
the functional form of the decay depends strongly upon the location of the initial quantum state with respect 
to the phase space of the underlying, unperturbed and perturbed, classical systems. In addition, the decay 
is sensitive to certain properties of the Hamiltonian perturbation. For instance, the Loschmidt echo decays 
differently depending on whether the perturbation is global or local in phase space, whether it vanishes under 
averaging over time, or whether it preserves or destroys integrability of the underlying classical system. A 
number of different regimes, such as the exponential decay and power-law decays with varies decay exponents, 
have been observed in numerical simulations. Interesting phenomena of quantum revivals, closely related to 
(quasi-)periodicity of the underlying classical motion, and temporary quantum freeze of the Loschmidt echo 
have been reported in various studies. However, a comprehensive classification of all decay regimes of the 
Loschmidt echo in regular systems still remains a challenge in the field. 

2.3.3 Numerical observation of various decay regimes 

Refs. |Jacq 01| , ICucc 02al ICucc 02bl TCucc 041 IGous 071 lGbus 081 lAres 091 iGarc lib] 

Numerical simulations have been of foremost importance for establishing the results summarized in sections 
12.3.11 and 12.3.21 The main body of numerical work has been devoted to one-particle systems. Various chaotic 
dynamical systems, globally perturbed, have been studied numerically; among them the Lorentz gas, two- 
dimensional hard-wall and soft-wall billiards. As a prominent example, Fig. [7] shows the crossover from the 
perturbative Fermi-golden-rule regime to the Lyapunov regime of the decay rate of the Loschmidt echo as a 
function of the perturbation strength in the Lorentz gas. The predictions for local perturbations, described in 
Sec. 12.3.11 were also tested in a chaotic billiard system. In this case, a deformation of a small region of the 
boundary was considered as the perturbation (see Fig. HJ). 

Other enlightening numerical studies were carried in quantum maps. These systems possess all the essential 
ingredients of the chaotic dynamics and are, at the same time, extremely simple from a numerical point of 
view, both at the classical and at the quantum levels. The predictions for locally perturbed chaotic systems 
were also observed in the paradigmatic cat map. Figure shows the decay rate of the Loschmidt echo as a 
function of the perturbation strength k in the cat map under the action of a local perturbation. When the 
map was perturbed in all its phase space, three decay regimes (instead of two) were observed depending on the 
perturbation strength: (i) the Fermi-golden-rule regime for weak perturbations, (ii) a regime of an oscillatory 
dependence of the decay rate on k at intermediate perturbation strengths, and (iii) the Lyapunov regime for 
strong perturbations. The complete understanding of the crossover between the last two regimes remains an 
open problem. 

2.3.4 Beyond the "standard" picture: fluctuations, disorder, many-body systems 
Refs. jAdam 031 ISilv 031 iPeti 051 IManf 061 |Quan~06| IFozzToTl IManf 081 |Zang| 

Decay regimes of the average Loschmidt echo, discussed above, are the ones most commonly observed, and 
generally regarded as "standard" , in low-dimensional quantum systems. However, when analyzing the Loschmidt 
echo beyond a simple averaging over the initial state or perturbation, or in a more complex setup, one often 
encounters departures from the standard picture. 
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Figure 7: Decay rate of the Loschmidt echo as a function of the perturbation strength k for the three concen- 
trations of impurities in a Lorentz gas. The rates F are normalized to the group velocity of the initial wave 
packet. The solid lines are the corresponding classical Fyapunov exponents and the dashed lines are fits to 
the quadratic behavior. Inset: a log-log scale of the same data showing the quadratic increase of T for small 
perturbations. Adapted from Ref. [Cucc 04j . Copyright (2004), American Physical Society. 

It was fist pointed out by Silvestrov and coworkers that, for sufficiently strong perturbations and on time 
scales short compared to the Ehrenfest time, statistical fluctuations generally play an important role in the 
problem of the Foschmidt echo decay. In particular, the average Foschmidt echo, M(i), is typically dominated 
by rare fluctuations, characteristic of only a small fraction of the chaotic phase space. In the main part of the 
phase space, however, the unaveraged Foschmidt echo, M(t), decays much faster than M(t): this decay can 
be as fast as double-exponential, M(t) ~ exp(— constant x e 2At ). It is only after the Ehrenfest time that M{t) 
follows M(t) in the main part of the phase space. 

The variance of the Foschmidt echo, M 2 — M , was addressed by Petitjean and Jacquod within both a 
semiclassical and a Random Matrix Theory approach. The variance was shown to exhibit a rich nonmonotonous 
dependence on time, characterized by an algebraic growth at short times followed by an exponential decay at 
long times. 

In disordered systems, the Foschmidt echo may exhibit even richer decay than that in low-dimensional dy- 
namical systems, since the elastic mean-free path and the disorder correlation length enter as relevant scales. 
For long-range disorder (small-angle scattering) an asymptotic regime governed by the classical Fyapunov expo- 
nent emerges. A Josephson flux qubit operated at high energies leads to a chaotic dynamics, and the Fyapunov 
regime can be obtained when a Foschmidt setup is considered in such a system. 

The Foschmidt echo studies in many-body systems are comparatively less developped than those in one- 
body cases. The numerical calculation of the many-body Foschmidt echo is a highly demanding computational 
task, and the standard approximations are in difficulty for describing the small difference obtained after the 
forward and backward time-evolutions. For this reason, the most reliable results of Foschmidt echo in many- 
body systems are those of a one dimensional spin chain or in an Ising model with transverse field, where the 
influence of criticality in the fidelity decay has been established. Trapped Bose-Einstein condensates have been 
treatead as a many-body Foschmidt echo setup using mean-field approaches. 

2.4 Phase-space representation and classical fidelity 
Refs. [Bene 02l Bene 03I lEckh 03I iGarc 03l ICucc 04l rVebToi] 

An intuitively appealing representation of the Foschmidt echo in phase space is obtained by expressing Eq. (TTJ) 



where Wh x and Wh 2 are the Wigner functions resulting from the evolution of IFo(q, p), introduced in Eq. (J7]), 
under the action of the Hamiltonian operators Hi and H2 respectively. This approach is naturally connected 



as 




(12) 
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Figure 8: The Loschmidt echo decay in the desymmetrized diamond billiard for four different values of the cur- 
vature radius r of the arc boundary deformation. The solid straight line gives the trend of the exp(— 2jt) decay. 
Inset: Forward-time wave packet evolution in the unperturbed billiard, followed by the reversed-time evolution 
in the perturbed billiard. The arrow shows the momentum direction of the initial wave packet. The propa- 
gation time corresponds to approximately 10 collisions of the classical particle. Adapted from Ref. [Gous 07] , 
Copyright (2007), American Physical Society. 

to the dephasing representation of Sec. 12.2.11 and provides a particularly useful framework in the study of 
decoherence, as the Wigner function is a privileged tool to understand the connection between quantum and 
classical dynamics. 

An initial Gaussian wave-packet is associated with a Gaussian Wigner function that will develop in phase- 
space non-positive structures and be deformed under the evolution of H± and H2 . The sensitivity of the overlap 
(Tl2"|) to the non-positive parts and to the deformations of the Wigner function can be related with the different 
decay regimes discussed in Sec. 12.3.11 

The form ([12"]) of the quantum fidelity allows to define a classical fidelity by using in a classical problem 
the Liouville distributions Lh 1 and Lh 2 instead of Wh x and Wh 2 respectively. This definition follows from the 
logical representation of the Liouville distribution as the classical limit of the Wigner function. However, the 
definition by an overlap of distributions does not convey the sense of classical reversibility: a trajectory evolving 
forward in time with Hi and backwards with Hi would give a considerable contribution to the overlap not only 
if it ends up close to the initial point, but whenever the final point is in the neighborhood of those of the initial 
distribution. The labeling of trajectories and particles, characteristic of classical mechanics, is then not taken 
into account in this definition of classical fidelity. 

The correspondence principle dictates that the quantum fidelity follows its classical counterpart up to the 
Ehrenfest time. However, such a correspondence does not imply that, in the intermediate-time asymptotic decay 
of globally perturbed chaotic systems, the Lyapunov regime is only present up to this characteristic time scale. 
In fact, the Lyapunov regime is not simply an effect of the classical-quantum correspondence. For instance, in 
a two-dimensional billiard, the saturation time is (see Sec. 12.3. ip that would signal the end of the asymptotic 
decay (and so the end of the Lyapunov decay if one is in the appropriate regime) can be much larger than the 
Ehrenfest time. 

3 Experiments 

3.1 Nuclear Magnetic Resonance 

Ref. [Slic 90] (see particularly Sec. 8.8 and App. E) 

Nuclear magnetic resonance has been the main tool to study echoes generated by different time-reversal pro- 
cedures since the fifties. After the time-reversal of individual spin precession implemented by Erwin Hahn in 
the spin echo, one had to wait until the seventies, when Rhim and Kessemeir achieved the reversion of the 
spin-spin interactions. The resulting reversal of the macroscopic polarization is known as Magic Echo and 
was thoroughly studied by Rhim, Pines and Waugh. In the nineties Ernst et al. implemented a strategy to 
address a localized spin excitation and Levstein et al. realized that such study was optimal for a theoretical 
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Figure 9: Decay T of the Loschmidt echo of a cat map perturbed locally with a shear in momentum as a function 
of the scaled perturbation strength n. The symbols correspond to the following fraction of the phase space that 
is perturbed: V (40%), A (30%), + (20%) and * (10%). The width 77 of the local density of states is also 
plotted with solid lines. This width gives the spread of the states of the unperturbed system in the bases of the 
perturbed one It is noticeable that, the smaller is the perturbed region, the more similar to rj V is. Adapted 
from Ref. |Ares 09] . Copyright (2009), American Physical Society. 



description. They used crystal structures with abundant of interacting 1 H spins, where the dipolar coupling 
produces the "diffusion" of the initial polarization. The time-reversal procedure is then able to refocus the 
excitation generating a Polarization Echo. These last experiments have a direct connection with the concepts 
discussed in this review. 



3.1.1 Basic concepts on spin dynamics 

Refs. [Past 951 iMadi 971 iLevs 981 iDanToil [Zang| 

A typical experiment starts with a sample which constitutes a network with about 10 23 interacting spins in 
thermal equilibrium in presence of an external magnetic field Bq. This many-spin state is denoted as |W cq ). 
In this state every spin has an almost equal probability of being up or down since the thermal energy, ksT, is 
much higher than any other relevant energy scale of the system. At time t = a sequence of radiofrequency 
pulses ensures that spin at site 0th is oriented along z direction. This is represented by the action of the spin 
operator Sq on \ip e q)- In a system with m + 1 spins, we represent this initially excited state as 



l*°> - |/T ,^!!| q l - E ^ Ito) a IA-) , (13) 
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where the denominator | (^ eq | 5(7 S^j ^eq) | ensures that the initially excited state has a proper normalization 
and 4> r is a random phase that describe a mixture of states of the form 

|/3r) = <8> \si) <8 |s 3 ) <8 ... <8 \s m ) with \s k ) e {|t) , II)} • (14) 

The preparation of this state involves the use of a 13 C nucleus as a "spy" to inject and detect polarization 
at the directly bonded 1 H spin (0th). Then, the system evolves under mutual many-body interaction described 
by Hi for a period t%. In the polarization echo experiment, this is described by effective dipolar interaction 
Hamiltonian Hi, truncated by the Zeeman field, 



Polarization £- ' 

Echo Experiment % A 



(15) 



Hi contains flip-flop or XY operators of the form S^~S~ + S~Sj~, as well as Ising terms of the form S*Sj. 
The dipolar interaction dij constant decay with the third power of the distance between sites i and j. Hi 
produces the spread of the initially localized excitation. Since total polarization is conserved under Hi , such 
process is commonly known as "spin diffusion" . Then, a new pulse sequence rotates all the spins 90 degrees and 
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the irradiation with r.f. field produces a further truncation of the dipolar interaction along the rotating field. 
The resulting effective Hamiltonian — H2 = — {Hi + E) acts for another period t%. Observe that the reversal 
is not perfect but E is a small term which is due to truncations required to form the effective Hamiltonian. 
This constitutes the "backwards" evolution period after which the local polarization at site Oth is measured. 
We notice that E acts as a self-energy operator that may account for interactions within the same Hilbert 
space, and thus is just an Hermitian effective potential. However, it also may describe the interaction with 
a different subsystem that constitutes the "environment". In the last case E = A — ir would also have an 
imaginary (non-Hermitian) component. In this last case, evolution of observables should be described by two 
coupled Lindblad or Keldysh equations for the density matrix requiring a self-consistent evaluation of dynamics 
of the system and the environment. This last situation would make the analysis less straightforward and might 
hinder the essential physics. Thus, we first focus on perturbations that ensure unitary evolution. An example 
discussed by Zangara et al. of this situation is a second frozen spin chain interacting with the first through an 
Ising interaction. In these cases, a polarization echo is formed when t% = t\, i.e. after a total evolution time 
2t = £2 + t\- The normalized amplitude is given by the spin autocorrelation function that accounts for the 
observed local polarization, 

(* cq |g -(2t) S + |*c 
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where 5^~(2t) = ^Hit/^Q-^t/hg- e iH 2 t/h e ~iH 1 t/h j g local spin operator, expressed in the Heisenberg rep- 
resentation respect to the acting Hamiltonian H(t) = H\0{t — t) — H^Qir — t). An equivalent definition of 
Mpe(£) in terms of S§(2t) Sq was also proved useful. 

The connection of this experimental observable with the concepts discussed so far are already hinted but 
not yet clarified by Eq. (1161) . Thus, we consider that the m+1 spins are arranged in a linear chain or in an odd 
sized ring where spin-spin interaction is restricted to XY terms acting on nearest-neighbors. In this case the 
Wigner- Jordan spin-fermion mapping allows to describe the spins system as a gas of non-interacting fermions 
in a 1-d lattice. In that problem, the initial excitation propagates as a single spinless fermion. In turn, the 
observed local polarization can be written back in terms of the evolution of an initially localized "spin wave" : 



\ipo) = Ito) ® 111) ® 14-2) ® lis) ® ■■• ® \U 
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where ||j n ) (|4-n)) means that the n-th site is occupied (unoccupied). Thus, it is clear that the excitation at 
site 0-th can spread along the 1-d chain through the nearest neighbor flip-flop interaction. This correspondence 
between many-body dynamics and spin wave behavior was worked out in detail by Danieli et al. This procedure 
was originally employed by Pastawski et al. in the context of time-reversal experiments to predict the presence 
of Poincare recurrences (Mesoscopic Echoes) which were clearly observed by Madi et al. at the laboratory of 
Richard Ernst in Zurich. In this condition, the polarization detected after the time-reversal procedure results: 
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Notice that here \tpo) denotes a single particle wave function, there are no other operators than the propagators 
and the modulus square makes explicit the positivity of Mpe(£). These features that were not obvious in 
Eq. (|16p and make it to agree with the definition of the Loschmidt echo of Eq. |T|) . 

Many other effective many-body Halmiltonians can be experimentally reversed to obtain the revival of 
different observables that reduce to a detectable polarization. Thus it is now a common practice to losely call 
Loschmidt Echo to the polarization resulting from any of these possible spin dynamics reversal procedures, 
regardless the applicability of the single particle correspondence, with an explicit mention of the observable 
used. 



3.1.2 Time-reversal of the dipolar many-spin interaction 
Refs. |Levs 981 [Usaj 98| , IPast 001 [Zure~07] 

The experiments have been carried out in different molecular crystals of the cyclopentadienyl (C5H5) family. 
There, 1 H nuclei are arranged in five fold rings with strong intra-ring dipolar interactions which nevertheless are 
not constrained to the molecule but also yield appreciable intermolecular couplings. Polarization can be injected 
on those 1 H which arc close to a "spy" 13 C nucleus that acts as a local probe used to inject and eventually detect 
local polarization as represented in Fig. 1101 Spreading of the initally localized spin excitation proceeds through 
a dipolar interaction which is truncated either in the basis of external magnetic field (Laboratory frame) or in 
the rotating frame of a radio frequency field. This choice is crucial in order to achieve the change of the sign and 
the eventual scale down of the effective Hamiltonian. Thus a spin dynamics with Hi in the Laboratory frame is 
reversed by a dynamics of Hi in the rotating frame or vice versa. The role of environment might be played by 
paramagnetic Co(II) or quadrupolar Mn nuclei that replace the Fe and thus they are absent for pure ferrocene 
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Figure 10: Ferrocene crystal showing the 1 H (small dots), C atoms (medium-sized dots) and Fe (large dots). A 
rare 13 C nucleus (usually called a "spy" nucleus), acting as a local probe, is shown by the arrow. Polarization 
evolves in the H lattice. Fe atoms at the center are not magnetically active, but when replaced by Co atoms 
relaxation of polarization is achieved. Adapted from Ref. [Levs 98] . Copyright (1998), American Institute of 
Physics. 

crystals. Besides, in pure ferrocene truncation in the rotating frame give rise to non-inverted non-secular terms 
which, being supressed by the corresponding Zeeman energy, are at most of the order of a few percent of the 
matrix elements of Hi- This would yield an Hcrmitian £. The strength of these terms is inversely proportional 
to B\, the strength of the r.f. pulse, and thus one has the possibility to experimentally reduce its importance 
by increasing the r.f. power. 

The results showing the buildup of the Loschmidt echoes, under H^ after different periods of evolution 
with Hi, are presented in Fig. 1111 The universal features of the Loschmidt echo appear when one studies 
the maximum recovered polarization (Loschmidt Echo intensity) as a function of t, the evolution time before 
reversal. In Ferrocene crystals, once one substrates the background noise it is clear that the Loschmidt echo 
follows a clear Gaussian law, 



for over two orders of magnitude and serves as a definition of the characteristic time T3. The same decay 
law is observed when, using a specially tailored pulse sequence, the Hamiltonians are scaled down by a factor 
n = 1,2,8 and 16 respectively, i.e. the characteristic decay time is increased by the facton n. The scaling of 
1 /T3 with the Hamiltonian strength is also observed when Hi and Hi are scaled down by using different crystal 
orientations. 

In contrast to Ferrocene, in Cobaltoncene crystals, one starts with a Gaussian decay, and as the dipolar 
Hamiltonian is scaled down, an exponential decay of the Loschmidt echo is developed and became stable with 
a characteristic time tse- Thus, the overall behavior of the recovered local polarization is 



One may assign this asymptotic exponential decay rate 1/tse to a Fermi Golden Rule decoherence rate induced 
by the perturbation induced by the paramagnetic nature of the Co (II) nuclei which acts as an environment. 
Further experiments by Pastawski et al. in crystals free of magnetic impurities confirmed the Gaussian decay 
with a rate I/T3 that depends only weakly on the r.f. power and, when truncation terms becomes too small, 
saturates to an intrinsic value that scales with the strength of the reverted Hamiltonian. 




(19) 
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Figure 11: Formation and fading away of the Loschmidt echo with reversal time ti elapsed after different 
periods t\ of forward evolution. Notice that the tails of the first Loschmidt Echo reflect the forward dynamics. 
At longer times there is a first revival asociated with the neighbors and a small revival associated with the 
Mesoscopic Echo resulting from the dynamics within the ring. Superimposed with the Loschmidt echo there are 
high frequency oscillations. These are coherent interferences resulting from an incomplete transfer between the 
13 C, used as a local probe, and the nearby 1 H nuclei. Adapted from Ref. |Levs 98j . Copyright (1998), American 
Institute of Physics. 

In summary, these experiments hinted that there is an intrinsic decoherence rate which is fixed by the 
inverted Hamiltonian. Thus, other than the difference between the Gaussian and exponential decay, one might 
say that I/T3 plays a role similar to the Lyapunov exponent in the reversal of chaotic one body systems. This 
surprising finding suggested that, even in absence of any important perturbation, in the experimental limit of 
m — > 00, even very small residual terms, in presence of the dynamical instability of a very complex many-body 
dynamics, are efficient enough to set the Loschmidt Echo decay into a perturbation independent decay regime. 
This plays the role of an "intrinsic decoherence" with a time scale determined by the inverted Hamiltonian. 

It was precisely the hypothesis of an "intrinsic decoherence rate" what triggered the theoretical analysis of 
time-reversal in chaotic systems. Indeed, while dynamics of one body systems in ID, described by Eq. (TT51) . 
can not be chaotic, one deems that disorder and the more complex many-body dynamics present in higher 
dimensional systems, have mixing properties which makes them assimilable to chaotic systems. This observation 
led to G. Usaj et al. to propose that quantum chaos contains the underlying physics of a time-reversal experiment 
in a many-body systems. It was also argued that Mpe(^) constitutes an entropy measure. These ideas finally 
boiled down into the model proposed by Jalabert and Pastawski. While Gaussian decays appears quite naturally 
as a consequence of the large number statistics of many-spin states that are progressively incorporated into the 
dynamics, the perturbation independent decay have not found yet a straightforward explanation in this context. 

3.1.3 Further time-reversal Experiments 

Refs. jLevs 04L [Sane 07L iRufe 09L ISanclM IRama llj 

While time-reversal of different interactions is almost an unavoidable tool in experimental NMR techniques, 
there are not so many studies devoted to grasp at the origins of the (in) efficiency of such procedures. In 
particular, time-reversal was implemented using different procedures and systems that range from 3c? and quasi- 
ld crystals to molecules in oriented liquid crystalline phases. In these cases, also different initial states and 
effective Hamiltonians were studied. 

Of particular interest are the double quantum (DQ) Hamiltonians, 



The notation d for the interaction strength recalls at it is an effective Hamiltonians built up from dipole-dipole 
interaction by suitable pulse sequences whose effect is described with the help of the average Hamiltonian theory. 
Since Hoq does no commute with the polarization described by Y^^S?, it gives a more mixing as dynamics 
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Figure 12: Attenuation of the polarization echo in the ferrocene single crystal as a function of t for n — 1 (REPE 
sequence) . The line represents a Gaussian fitting. The inset shows the PE attenuation for progressively reduced 
dipolar dynamics, n = 1, 2, 8, and 16. No asymptotic regime is reached within the experimental timescale. The 
solid lines are Gaussian fittings. Adapted from Ref. |Usaj 98] . 



explores different subspaces of the Hilbert space inducing the multiple quantum coherences. Besides Hjjq, 
there are high order terms that act as an environment whose instantaneous decoherence could be described by 
a Fermi Golden Rule. One experimental observation is that the larger the number of subspaces correlated by 
the interaction (coherences of high order) the faster the decoherence. 

One of the interesting properties of Id systems, is that the dynamics induced by a Hdq has a precise 
correspondence with an XY chain. Thus, the Loschmidt Echo is given by Eq. (|18p . In real systems, such 
as Hydroxyapatite (HAp), inter-chain interactions are significative, quantum Zeno effect induced by a strong 
dynamics along the chain prevents the development of coherences beyond second order and Eq. (|18p remains 
valid. This explains why the Loschmidt echo observed by Rufeil-Fiori et al. for the double quantum Hamiltonian 
shows a Loschmidt echo with an exponential decay described by the Fermi golden rule. 

This decay contrasts with results for an Adamantane crystal, a highly connected 3d system. There, a Fermi- 
function like decay is indicative of the two time scales. Each molecule has 16 spins which do not have direct 
interactions, thus remain independent for short times until they interact through neighbor molecules. In this 
case decoherence remains weak while intermolecular correlation builds up. Once neighboring molecules become 
fully coupled the density of directly connected states becomes very high and the Fermi-golden-rule controlled 
exponential decay takes over. 

Spin dynamics in various liquid crystals was studied in series of papers whose aim was to reduce the number 
of interacting spins to those within each molecule. However, the experiences have shown that a number of 
residual interaction remain significative and this strongly compromises the effectiveness of the Loschmidt echo 
sequences. 

3.2 Microwave billiards 

Refs. |Scha 05al IScha 05bl iHohm 081 |Body 09|, IKobe 101 [Kobe llj 



Microwave-frequency electromagnetic waves in quasi-two-dimensional cavities arc effectively governed by the 
Helmholtz equation, which is equivalent to the stationary Schrodinger equation. This equivalence allows one to 
address the Loschmidt echo in laboratory experiments with microwave billiards. 

Quantities directly measured in microwave experiments are frequency-dependent scattering matrix elements, 
S a b{v) and S' ab (v), corresponding to the unperturbed and perturbed chaotic billiard systems, respectively. Here, 
the indices a and b refer to the antennae (or scattering channels) involved in the experiment. The change from 
the frequency domain to the time domain is achieved by the Fourier transforms, S a b(t) = J dv e~ 2mvt S a b(v) and 
S' ab (t) — j dv e~ 27rlvt S' ab (v) , performed over an appropriate frequency window. The sensitivity of a scattering 
process to system perturbations are then quantified by the scattering fidelity amplitude 

f ab (t) = , , (22) 

(\S ab (t)\ 2 )(\Sa b (t)\ 2 ) 
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Figure 13: Attenuation of the polarization echo in the cobaltocene sample as a function of t for n = 1, 2, 5, 8, 
16. The experimental data show a clear crossover between a dominant Gaussian attenuation to an exponential 
one. The solid lines represent fittings of the whole set of data to Eq. (f2"U|) . Adapted from Ref. |Usaj 98| . 

where the asterisk stands for the complex conjugation, and the angular brackets denote an ensemble average over 
different realizations of the experiment (e.g., different positions of scatterers, antennae, etc.). The scattering 
fidelity itself is a real- valued quantity defined as F a b{t) = |/ a &(^)| 2 - The numerator in the right-hand-side of 
Eq. (|22[) quantifies correlations between scattering matrix elements of the unperturbed and perturbed system. 
The decay of the numerator however is dominated by the decay of the autocorrelations, so that the denominator 
is introduced to compensate for the autocorrelation contributions to the correlation decay. 

In chaotic systems and in the case of a weak coupling of the measuring antennae to the system, the scattering 
fidelity is known to approach the Loschmidt echo for a random initial state. This makes microwave billiards well 
suited for experimental studies in the field. In particular, microwave studies provided compelling experimental 
evidence for the validity of some of the known decay regimes of the Loschmidt echo for global and local 
Hamiltonian perturbations. 

3.3 Elastic waves 
Refs. jLobk 03L TGori 06j 

Sound waves that travel through a elastic medium are multiple scattered by its inhomogeneities generating a 
slowly decaying wave diffusion. If no change occurs in the medium over time, the usually called coda waves are 
highly repeatable, that is, for identical excitations the waveforms are indistinguishable. While if the medium 
is perturbed, the change in the multiple scattered waves will result in an observable change in the coda waves. 
Lobkis and Weaver have measured the sensitivity to temperature changes of elastic coda waves in aluminum 
alloy blocks. They used the cross correlation function between two signals obtained at different temperatures 
Ti and Ti 

X(e) = JMStAOStMI+s)) (23) 

y/fdtS^t) JdtS^(t(l + s)) 

to evaluate the distortion that is defined as D(t) = — ln(X max ), where the time dependence is given by the 
"age" of the signal. If formulated as a scattering process, it is shown that D(t) = — In [/(£)], where f(t) is 
the scattering fidelity that was introduced to analyze the fidelity decay in microwave billiards. For sufficiently 
chaotic dynamics in systems weakly coupled to decay channels, the scattering fidelity approaches the standard 
fidelity amplitude (the Loschmidt echo is the absolute square of the fidelity amplitude). 

The results obtained in acoustic signals traveling in an aluminum blocks were understood using the random 
matrix predictions for the standard fidelity (see Fig. [T7]). A surprising and unexpected finding is that the 
scattering fidelity decay for blocks with chaotic and regular regular classical dynamics are explained by the 
same random matrix expressions. 

3.4 Cold atoms 

Refs. jAnde 03L fAnde 04L lAnde 06L IWimb 06L [Mart 08L IWu~09l iKrTvTfl lUlla 111 IDube 12) 
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Figure 14: Loschmidt Echo experiment based on Hdq and -Hdq, for hydroxyapatite (HAp) (quasi-lci crystal) 
and Adamantanc, a fully 3-d crystal. Here, the source of £ are very small high order terms in the effective 
Hamiltonian. 

The Loschmidt echo, or rather a quantity closely related to the Loschmidt echo, has been carefully studied 
in laboratory experiments with cold atoms trapped inside optical cavities. These experiments offer an atom- 
optics realization of the echo in two-dimensional quantum billiards with underlying chaotic or mixed classical 
dynamics. The basic idea underlying the experiments with atom-optics billiards is as follows. An effectively 
two-level atom, with the internal states denoted by |1) and |2), is initially prepared in a state |1)<S> where 
stands for the spatial component of the initial state. A laboratory observation of the echo involves hrst exposing 
the atom to a sequence of microwave-frequency electromagnetic pulses and then measuring the probability Pi 
of finding the atom in the internal state |2). 

In the protocol targeting the Loschmidt echo, | (ip\e lH2t / h e~ lHlt / h \ip) | , an atom is hrst irradiated with a 
so-called n/2 pulse, which changes the atomic state into an equiprobable superposition of |1) ® and |2) ® \ip). 
The pulse is practically instantaneous and introduces almost no change to the spatial component of the state. 
Then, the atom is left to evolve in an optical trap for a time t, during which components |1) and |2) of the atomic 
state propagate under Hamiltonians H\ and H% respectively. The difference between H\ and Hi originates from 
a difference in dipole interaction potentials exerted upon the states |1) and |2) by the optical trap. Finally, 
another 7r/2 pulse is applied to the atom and the probability Pi of finding the atom in the internal state |2) is 
measured; Pi turns out to be a quantity closely related to the Loschmidt echo. 

Realistic experiments however deal with thermal incoherent mixtures of initial states, rather than with pure 
states. As a result, the echo amplitude, due to an individual state of this thermal mixture, contributes to Pi 
with an effectively random phase, smearing out the echo signal. In order to overcome this difficulty, Andersen, 
Davidson, Griinzweig, and Kaplan addressed a new echo measure, |^| e l #2V 2 fte lHlt / 2ft e _lff2t / 2ft e _lHlt / 2fi |^>}| , 
closely related to the Loschmidt echo. In their experiments, "two-level" atoms, initially prepared in the state 
|1) <8> \ip), were successively exposed to (i) a ir/2 pulse, (ii) evolution for a time t/2, (iii) a ir pulse swapping the 
population of the internal states |1) and |2), (iv) another evolution for a time t/2, (v) a tt/2 pulse, and finally 
(vi) a measurement of Pi . Under this protocol, each state of the thermal mixture contributed to the probability 
Pi with the same phase. This allowed for a reliable measurement of the echo even for ensembles of more than 
a million of thermally populated states. 

More recently, atom interferometry has also been used to investigate several aspects of time-reversal in 
quantum kicked rotor systems. Wu, Tonyushkin, and Prentiss studied the dynamics of laser-cooled rubidium 
atoms subjected to periodically pulsed optical standing waves, as an atom-optics realization of the kicked rotor. 
Their experiments have demonstrated that quantum fidelity of a system, that is chaotic in the classical limit, 
can survive strong perturbations over long time without decay. In a similar setup, Ullah and Hoogerland have 
performed an experiment that demonstrated the possibility of using time-reversal evolution for cooling atomic 
matter waves. The problem of fidelity decay in kicked atom-optics systems has recently drawn considerable 
interest among theoreticians. 

3.5 Time-reversal mirrors 
Refs. jFink 991 ICalv 081 iTadd 09] 
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Figure 15: Correlation function C[Sii, S'n] oc (S^t) S'uft)) . The experimental results for the autocorrelation 
are shown in black, while the correlation of perturbed and unperturbed system are shown in gray/orange. The 
smooth solid curve corresponds to the theoretical autocorrelation function, and the dashed curve to the product 
of autocorrelation function and fidelity amplitude. The inset shows the billiard geometry used. Movable parts 
are marked with an arrow. Adapted from Ref. |Scha 05b] . Copyright (2005), American Physical Society. 

Another time-reversal attempt conceptually close to the Loschmidt echo is that of time-reversal mirrors, de- 
veloped in the last twenty years. Such a procedure has been successfully implemented in various setups where 
classical waves propagate through a complex media (going from acoustic to electromagnetic waves). In the 
time-reversal mirror protocol an initially localized pulse is recorded by a collection of receiver-emitter transduc- 
ers during a time interval where the wave suffers multiple scattering. The later re-emission of the time-reversed 
signal by each transducer during an interval of time equal to the recording one leads to the refocusing of the 
signal in the region of the original excitation. 

Two features of this protocol have lead to considerable surprise among the practitioners: 

(i) even one transducer is enough to obtain a good reproduction of the original signal 

(ii) the quality of the refocusing was improved by the complexity of the media yielding the multiple scattering 
of the waves. 

A semiclassical theory of time-reversal focusing can be built up in terms of propagators and classical trajec- 
tories. The two previous features can be acknowledged by the semiclassical theory. 

Time-reversal mirrors and Loschmidt echos differ since the former aims the refocalization of a wave which 
is localized in space and time, while the latter attempt to time-reverse a quantum state. A common aspect 
between both protocols is the fact that the Hamiltonian for the forward and backward evolutions can differ due 
to modifications of the environment during the process. 

Time-reversal mirrors are not only conceptually important, but also have very important technological 
applications as for example brain therapy, lithotripsy, nondestructive testing and telecommunications. Recently 
a sensor of perturbations was proposed and demonstrated combining the ideas of the Loschmidt echo and time 
reversal mirror in classical waves. 

4 Past, present, and future of the Loschmidt echo studies 

Since the early discussions between Boltzmann and Loschmidt on irreversibility and time-reversal, it has been 
clearly established that chaos is the source of the irreversibility in statistical mechanics. However, when con- 
sidering the sources of irreversibility in a quantum system with a few degrees of freedom Quantum Mechanics, 
which is the fundamental theory of the microscopic world, does not allow for chaotic behavior in the sense in 
which the latter appears in Classical Mechanics. That is, there is no exponential separation between states with 
nearby initial conditions since quantum evolution is unitary. Trying to understand the origin of irreversibility in 
quantum mechanics, Asher Peres proposed in 1984 as an alternative, to study the stability of quantum motion 
owing to perturbations in the Hamiltionian. In his seminal paper, Peres considered M(t), later called Loschmidt 
echo, as a measure of sensibility and reversibility of quantum evolutions. He aimed to distinguish classically 
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Figure 16: Decay rate of the Loschmidt echo in a microwave billiard as a function of the displacement h of 
a small flat part of the billiard boundary. The solid (dashed) curve corresponds to the decay rate predicted 
by the semiclassical theory taking (not taking) into account classical trajectories multiply scattered within the 
boundary deformation region. Adapted from Ref. [Kobe llj . Copyright (2011), American Physical Society. 

chaotic or integrable dynamics according to the speed at which the fidelity decays. Peres reached the conclusion 
that the long-time behavior of the fidelity (or saturation, following the terminology of Sec. 12.3.11) in classically 
chaotic systems is characterized by smaller values and smaller fluctuations as compared to the case of regular 
dynamics. In view of numerous possibilities of different behaviors that are allowed in regular systems, such a 
conclusion does not always hold. 

It is interesting to notice that Peres' article appeared almost at the same time than two other seminal works 
studying the relevance of the underlying classical dynamics for quantum properties. These are the random- 
matrix description of the statistical properties of classically chaotic systems by Oriol Bohigas and collaborators; 
and the understanding of the spectral rigidity from a semiclassical analysis of periodic orbits proposed by 
Michael Berry. The field of Quantum Chaos developed building upon these founding ideas and most of the 
subsequent works were concerned with the spectral properties of classically chaotic systems. Comparatively 
very little work was done along the Peres' proposal in the decade following it. 

Motivated by the puzzles posed by echo experiments in extended dipolar coupled nuclear spin systems, 
Rodolfo Jalabert and Horacio Pastawski studied in 2001 the behavior of the Loschmidt echo for classically chaotic 
systems. They found that, depending on the perturbation strength, the decay of the Loschmidt echo exhibits 
mainly two different behaviors. For weak perturbations, the decay is exponential with the rate that depends on 
the perturbation strength and that is given by the width of the local density of states (usually called the Fermi- 
golden-rule regime). For stronger perturbations, however, there is a crossover to a perturbation-independent 
regime, characterized by an exponential decay with the rate equal to the average Lyapunov exponent of the 
underlying classical system. 

The connection of the quantum Loschmidt echo with the classical chaos generated a great activity in the field. 
Researches from different fields, such as quantum chaos, solid state physics, acoustics, and cold atom physics, 
have made important contributions towards understanding various aspects of the Loschmidt echo. During the 
first years of the last decade the studies were mainly focused on one-body aspects of the problem and on the 
connection between the Loschmidt echo and decoherence; many interesting experiments were also performed. 

In the last years, the interest has primarily been focused on various many-body aspects of the Loschmidt 
echo. For example, there are studies that consider the decay of the Loschmidt echo as a signature of a quantum 
phase transition, or that concentrate on the relation between the Loschmidt echo and the statistics of the work 
done by a quantum critical system when a control parameter is quenched. In view of recent improvements in 
experimental many-body system techniques, the issue of the Loschmidt echo decay has become more concrete. 

Some of the more important open questions in the field are: 

• The Lyapunov (perturbation-independent) decay regime has been one of the most influential break- 
throughs in the theory of the Loschmidt echo. However, an experimental observation of the regime 
in simple well-controlled systems is still lacking. 
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Figure 17: The distortion as a function of time for coda waves that travel in an aluminum alloy block at different 
temperatures. The thin lines correspond to measurements in different frequency ranges. The thick lines show the 
corresponding theoretical curves obtained using Random matrix theory. Inset: Laboratory configuration. An 
aluminum specimen, typically with nonparallel faces and a defocusing cylindrical hole, is cooled in a vacuum as 
temperature and ultrasonic response are monitored. A reed relay isolates the response from the pulser. Adapted 
from Refs. [Lobk 03] and [Gori 06] . Copyright (2003, 2006), American Physical Society. 

• The first experimental measure of the Loschmidt echo was done in a many-body spin system using NMR 
techniques. However, theoretically little is know about the behavior of the Loschmidt echo in many-body 
systems. Very interesting experiments are now being carried on using NMR. 

• The semiclassical theory of the Loschmidt echo has proved to be a powerful tool to understand many 
aspects of its behavior. But this approach has some well-recognized difficulties, such as the root search 
problem or the exponentially growing number of classical orbits needed for the semiclassical expansion. 
Recently, a simple semiclasical dephasing representation of the Loschmidt echo was proposed that does 
not suffer of the usual problems of the semiclassical theories. The range of validity is however unknown. 

• There is a vast amount of work characterizing the decay regimes of the Loschmidt echo as universal. 
For chaotic systems perturbed with global or local perturbations the behavior of the decay rate of the 
Loschmidt echo is well established (see Sec. 12.31) . However, recent works in quantum maps that are 
perturbed in all phase space, have reported non-universal oscillatory behavior in the decay rate of the 
Loschmidt echo as a function of the perturbation strength. Deviations from the perturbation indepen- 
dence are found usually in the form of oscillations around the Lyapunov exponent. Moreover, there are 
cases where deviations are considerably large rendering the Lyapunov regime non-existing. The fully un- 
derstanding of this non-universal behavior, that was also observed in a Josephson flux quit, is still lacking 
[Garc llbllFozz 07] . 

• The semiclassical theory of the Lyapunov decay of the Loschmidt echo is based on highly localized initial 
states and the diagonal approximation of Eq. (3) for M(t). But recent results using the semiclassical 
dephasing representation have shown a Lyapunov regime in the mean value of the fidelity amplitude 
(m(t)} [Garc llallZamb 11] , 

• Systems with regular or fully chaotic dynamics are rather exceptional in nature. A generic system has 
mixed phase space consisting of integrable islands immersed in a chaotic sea. Few general results are 
known about the Loschmidt echo in this scenario. 
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